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The high energy evolution equations that describe the evolution of hadronic amplitudes with 
energy are derived assuming eikonal interaction of the evolved hadronic wave function with the 
target. In this note we remark that this derivation allows a different interpretation, whereby the 
hadronic wave function is not evolved, but instead the evolution acts on the S - matrix operator. 
In this approach, analogous to the Heisenberg picture of Quantum mechanics, the scattering is not 
eikonal and additional boost provides for radiation of more gluons in the final state. 

This short note continues the program of investigating high energy evolution of hadronic scattering observables[l], 
[2], [3], [4], [5], [6], [7], [8], [9], [10]. In the past several years the activity in this area concentrated on attempt- 
ing to include the effects of Pomeron loops in the evolution with many contributions by several groups 
[11], [12], [13], [14], [15], [16], [17], [18]. In particular a recent paper [19] gives a derivation of the evolution of a hadronic 
light cone wave function valid for all parametrically interesting values of the valence color charge density, from the 
weak charge density (KLWMIJ) limit to the large charge density (JIMWLK) limit. How to use this information 
to derive the evolution of scattering amplitude is still an open question. The physical picture of the evolution as 
described for example in [20], is that the increase in energy is used to boost the projectile wave function. The boost 
leads to increase of the number of gluons in the projectile wave function. These gluons then scatter eikonally when 
propagating through the color field of the target. 

The statement that the high energy scattering in non abelian theories is eikonal frequently causes eyebrows to lift, 
as it is well known from diagrammatic calculations that at high energy logarithmic contributions to the S matrix 
come from final states with extra gluons. The confusion here stems from the fact that even though each parton 
scatters eikonally, the whole projectile state of course does not and so inelastic final states certainly contribute to the 
scattering. The evolved wave function contains more partons. Those partons decohere while propagating through the 
target field, even though individual partons do not emit gluons. Due to this decoherence the majority of final states 
are inelastic. We point out in this note that the derivations of [9] and [14] can be also given another interpretation. 
This is akin to the Heisenberg picture of quantum mechanics. Here it is not the projectile state that is evolved, but 
rather the ^-matrix operator. The S matrix thus has explicitly noneikonal contributions which are responsible for 
emission of extra gluons into the final state. 

Consider calculation of the expectation value of any observable in the wave function of the incoming probe \v) . We 
will be interested only in observables which depend on the color charge density operator j a (x) and not on any other 
property of the wave function 

= (v\0[j]\v) (1) 
An example of such an observable is the eikonal S - matrix in the external field a 

S = e X p{ijj a (x)a a (x)\ (2) 

In the approach of [9] the calculation of the scattering amplitude for a hadron-hadron scattering is given by 

S = (S) a (3) 

where the weight for the averaging over the fields a is furnished by the target wave function. 

The projectile wave function at certain rapidity contains partons with longitudinal momenta larger than a certain 
cutoff A. As discussed many times in the literature, the standard way to evolve the observable with rapidity is to 
boost the projectile wave function. Under boost, the longitudinal momenta of the existing partons are scaled by the 
Lorentz 7 factor, while the newly opened phase space is filled by soft gluons. 

\v)^n Y [j,a\a]\v) (4) 

The operator Vty creates the soft gluons in the projectile wave function. It is a Bogoliubov type operator and has 
been calculated in [19] for arbitrary charge density j. In the low density limit, that is when j = O(g) the form of this 
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operator has been found in [14]. In this case it reduces to the coherent operator 



fV-CV=exp{2i J d?xb«{x) J ^^-L m [al«(k+, X )+a«(k + ,x)] 



(5) 



where the integration is over the phase space opened by boost to rapidity Y, and the "classical field" b is determined 
by the color charge density through 



and the soft gluon creation and annihilation operators satisfy 

x), a{ f (p+ y)] = (2tt) 6 ab & S(k+ - p+) S^(x - y) . 
The observable O calculated in the state boosted to rapidity Y is given by 

Y = (v\n Y 0\j a (x)+g [ ^-a\k+,x)T a a{k + ,x)]n Y \v) 



(6) 



(7) 
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with T£ 



if abc . The charge density operator j is shifted here since the newly produced gluons also carry color. 
Since we are interested in the dependence of observables on rapidity, the rapidity in this approach plays the role 
of time - the parameter of the evolution. The operator f2 then plays the role of the evolution operator in rapidity. 
The state \v)y — Qy\v) is analogous to the time dependent state in the Schroedinger representation of Quantum 
mechanics. Thus the picture of the evolution of the scattering amplitude for example, is that the projectile state is 
evolved while the scattering matrix operator always stays eikonal. There are more partons (gluons) in the evolved 
state, and the scattering matrix element changes with rapidity, even though the partons scatter eikonally and do not 
emit any extra gluons during the propagation through the target. 

One can give an alternative interpretation to eq.(8). In particular we can think of the operator Qy as acting on the 
observable O rather than on the state \v). In this interpretation it is the operator O that is evolved to higher rapidity 
rather than the state \v). This is analogous to the Hcisenberg picture of Quantum mechanics. Thus, for example 
if at low enough rapidity we take the operator O as eikonal scattering matrix eq.(2) its evolution with fl produces 
a new S matrix operator which is not eikonal anymore, but rather allows for emission of final state gluons. This is 
straightforward to see using the explicit form of the action of the operator on the fields which has been given in 
[19]. For simplicity instead of general ft we will use here its low density version C. Thus to order g 2 we have 



J a (x) 



■a){k+,x)T a a{k+,x) 



^-a^(k+,x)T a a(k+,x) 
-2igb»(x)lZA$(x, 0) - 2i±J d 2 y^^j b (x) P' V Ml . ())] 
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where we have defined 



dk+ 1 



-.[a^k + ,x) + a\ c (k + ,x)] 



(9) 



(10) 



Obviously the "rapidity evolved" color charge density contains terms linear and quadratic in the soft gluon creation 
operator. Thus the "rapidity evolved" S- matrix 



S Y = C f SC = 



exp 



i / C*j a (x)Ca a (x) 



(11) 



describes not only eikonal propagation of original partons, but also emission of soft gluons by these partons while 
they propagate through the target field. As is obvious from eq.(9) the amplitude of such emission is proportional to 
the charge density in the projectile, the field strength in the target and the rapidity interval Y over which the process 
has been evolved. 
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We note that the original derivation of the JIMWLK evolution equation that describes the high density limit was 
in fact performed by calculating the change of the color charge density due to the evolution. In this sense the original 
derivation uses the Heisenberg picture of the evolution. 

The derivation of the KLWMIJ equation on the other hand was given in the Schroedinger picture. In the rest of this 
note we will rederive this equation in the Heisenberg picture. The purpose of this exercise is twofold. Firstly, we want 
to put both derivation in the same framework. Secondly and more importantly we anticipate that the generalization 
of the high energy evolution to the general case (beyond JIMWLK and KLWMIJ limits) will be more conveniently 
done in the Heisenberg picture. The reason is that in [19] the action of the operator O on all the degrees of freedom 
of the theory has been calculated. Thus we know the explicit form of the evolved color charge density. On the other 
hand the explicit form of the operator £1 in terms of dynamical fields has not been found, and therefore its explicit 
action on the wave function \v) is not available. Thus we want to make our hands a little dirty and gain some more 
experience working with the Heisenberg representation. 

Any observable 0[j] can be expanded in Taylor series in powers of j. Thus it is enough to consider 



When evolved to small rapidity AY 



O n =r{ Xl )...j an {Xn) 



O n -» 0\ Y = C\j ai { Xl )...J a -{x n )C 



with 



/dk + 
—a\k+,x)T a a(k+,x) 



(12) 



(13) 



(14) 



Our aim is to represent the expectation value of the evolved operator O n in the valence state in terms of a 
" Hamiltonian" acting on the expectation value of the unevolved operator 



(v\Ol Y \v) = (v\O n \v) + AY(v\H 



O n \v) 



(15) 



The above equation is somewhat little cryptic, since j is a quantum operator and differentiation with respect to 
quantum operator has to be defined. To make its meaning precise we remind the reader that averaging over the 
valence wave function can be done in the path integral representation[14]. For this purpose one introduces an ordering 
variable t and considers classical functions j a {x, t). The average is then written as 

(v\j ai (x 1 )...j a »(x n )) = J DjW\j(x,t)]j^(x 1 ,t 1 )...j a "(x n ,t n ) (16) 

such that < t\ < t<i < ... < t n < 1. The values of the variable t are ordered in the same way as the position of the 
appropriate charge density operator in the operator product. This is necessary to account for noncommutativity of 
the quantum charge density operators. The weight function W is determined by the valence wave function. It is not 
arbitrary but has to satisfy particular constraints. This has been discussed at length in [17]. For our purposes it is only 
important to know that such a representation exists. The exact values of the time variable ti are also unimportant as 
long as the ordering is preserved. Given this the evolution eq.(15) takes the form of a functional differential equation 



W\j] 



(17) 



In eq.(17) j is a classical function which depends on the transverse coordinates x as well as the ordering coordinate t 
and functional differentiation is now well defined. 
To order g 2 we have 



{v\C*J ai (x 1 )...J a »(x n )C\v) = 

n—1 r „ 

J2 (v\3 a Hxi)-J am -Hx m -i) 2* / V a {x)A${x, 0),J a -(x m ) 
i L J x 



(18) 



m<k,m— 1 



J°- +1 (x m+1 )...J a ^(x k ^) x 



2i \ bi(x)A^x,0),J a "(x k ) 



(x k+1 )...j an {x n )\v) 



-if>ir^i)...r m -^ m -i)[2* f 61(^)^(^,0), [2* f 

m=l I Jx L Jy 



2i I bi(y)A b J (y,0),J a -(x r , 
y 



"- +1 (x m+ i 



...f-{x n )\v) 
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Here the charge density operators acting on the state \v) either to the right or to the left contain only valence part 
since the soft gluon contribution to charge density vanishes in the state \v) which contains no soft gluons. Let us 
hrst concentrate on the double sum term in eq.(18). First, we have to average over the soft gluon Hilbert space. 
This is quite straightforward. Only the term involving the annihilation operator a contributes in A(x, 0) in the first 
commutator, and only the term with contributes in the second commutator. The soft charge density operators 
between these two insertions pick up the contribution from the one gluon state created from \v) by the a> term, and 
thus are simply shifted by the one gluon charge. Once we have averaged over the soft gluons, the rest of the average 
can be represented as a functional integral of the form eq.(16). We can thus write 



n-l 



--Y I DjW\j] 

m<k,m—± \ 

f bl(x, t m - e)R ab {x, tm -e,t k - t)b\(x, t k - e) - f b\(x, t m + e)R ab (x, t m +e,t k - e)b[{x, t k - e) 

- f b l a {x,t m -e)R ab {x,t m -e,t k + e)b\{x,t k + e) + f VJx, t m + e)R ab (x, t m + e, t k + e)bl(x, t k + e) I x 

J X J X \ 

xj ai (x 1 ,t 1 )...j am - 1 (x m _ 1 ,t m _ 1 )j am {x m ,t m )..J ak {x k _ 1 ,t k _ 1 )j ak {x k (19) 
In this expression the time variables are strictly ordered t 1 < t 2 ... < t n . The operator 



R(x,t,t')=Vexp{g [ dr J T a 

{ Jt Sj a {X,T) 



(20) 



shifts the color charge density by the charge of a single gluon. The time coordinate of the field b is shifted by an 
infinitesimal amount e with respect to the coordinate of the nearest charge density so that the ordering of time 
coordinates follows the operator ordering of eq.(18). The differences between two pairs of terms combine into time 
derivatives if the summations are changed to time integrals so that we can write: 

J DjW\j]J dtj dt>^ J^(x,t)R ab (x,t,t/)b b i (x,t')r(x u t 1 )...j an (xn,t n ) (21) 

Although we have explicitly considered only the first term in eq.(18), it is straightforward to see that the second 
(double commutator) term reproduces the t = t' contribution in eq.(21). Integration over t and t' is straightforward 
with the result 



H 



r 6 ' 


i 







2 f bUx,0)R ab (x,0 7 l)b b (x,l)-K(x,0)bUx,0)~bt(x,l)bt(x,l) 

J x 



(22) 



This is precisely the KLWMIJ Hamiltonian first derived in [14]. 

This concludes the derivation of the KLWMIJ Hamiltonian in the Heisenberg representation of the evolution. We 
hope that this exercise will be useful for future derivation of the general evolution Hamiltonian including the Pomeron 
loop contributions using the results of [19]. 
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